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, Abstract. In self-gravitating stars, two dimensional or geophysical flows and in plasmas, long range 

1 ■ interactions imply a lack of additivity for the energy; as a consequence, the usual thermodynamic 
' limit is not appropriate. However, by contrast with many claims, the equilibrium statistical mechan- 
, ics of such systems is a well understood subject. In this proceeding, we explain briefly the classical 
■ approach to equilibrium and non equilibrium statistical mechanics for these systems, starting from 
I first principles. We emphasize recent and new results, mainly a classification of equilibrium phase 

a transitions, new unobserved equilibrium phase transition, and out of equilibrium phase transitions. 
^ ' We briefly discuss what we consider as challenges in this field. 

' O ' Keywords: Long range interactions. Inequivalence of ensembles. Kinetic theory. Out of equilib- 

Ch \ rium phase transitions. Two dimensional turbulence. Geophysical flows. 
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^ ■ 1. INTRODUCTION 

> • 

o : 

, In a large number of physical systems, any single particle feels a potential dominated 

^ I by interactions with far away particles: this is our definition of long range interactions. 

. ■ In a system with algebraic decay of the inter-particle potential V (r) ~ r", this occurs 

O ' when a is less than the dimension of the system (these interactions are sometimes called 

. "non-integrable"). Then the energy is not additive, as the interaction of any subpart of 

'. the system with the whole is not negligible with respect to the internal energy of this 

,>] given part. 

^ ■ Self gravitating stars, after the discovery of negative specific heat in have played 

a very important historical role, by emphasizing the peculiarities in the statistical me- 
chanics of systems with long range interactions. Besides astrophysical self gravitating 



systems [HjIEBSSIvLBSHIIII], the main physical examples of non-additive sys 
tems with long ran ge i nteractions are two-dimensional or geophysical fluid dynamics 
E Q H [iljSlivl] and a large class of plasma effective models JTsLM 0, 111. 



Spin systems 112211 and toy models with long range interactions [|23L have also 



been widely studied. The links between these different subjects have been emphasized 



recently Il23h . 



In these systems, the most prominent and interesting physical phenomenon is the self 
organization of the particles, or of the velocity field. This leads to coherent clouds of 
particles in plasma physics, to galaxies and globular clusters in astrophysics and to large 
scale jets and vortices in two dimensional or geophysical flows. Given the large number 
of particles or of degrees of freedom, it is tempting to adopt a statistical approach in 
order to describe these phenomena. The statistical description of such a self organization, 
both at the levels of equilibrium situations and relaxation towards equilibrium (kinetic 
theories), is a classical, long studied field. One of the aim of this proceeding is to insist 
on the vitality of this old subject and to stress new advances and remaining issues. By 
contrast, the out of equilibrium statistical mechanics of such phenomena is still in its 
infancy, and few studies have been devoted to it. We emphasize the importance of such 
studies for real applications, as most plasma and geophysical physical phenomena are 
out of equilibrium. We also describe some recent very suggestive results. 

Both equilibrium and out of equilibrium phase transitions play a key role in our under- 
standing of physics, because they separate regions of parameter space with qualitatively 
different behaviors. Very naturally, a large part of our studies will be devoted to phase 
transitions. We will especially stress the peculiar association of phase transitions with 
negative specific heat and statistical ensemble inequivalence in systems with long range 
interactions. We also insist on recently observed out of equilibrium phase transtions, in 
the context of two dimensional flows. Finally, we describe our personal guesses for what 
may be the challenges and interesting issues in the field of systems with long range inter- 
actions. We hope this could open new discussions, although we are conscious that such 
guesses are necessarily biased by personal prejudices. We actually sincerely hope that 
future researches will be much richer than what we describe. The article is organized 
in three main sections: equilibrium, relaxation to equilibrium and kinetic theories, non 
equilibrium stationary states. 

Equilibrium. Long range interacting systems are known to display peculiar ther- 
modynamic behaviors. As additivity is often seen as a cornerstone of usual statistical 
mechanics and thermodynamics, it is sometimes written in textbooks or articles that 
"statistical mechanics or thermodynamics do not apply to systems with long range inter- 
actions". In this paper, we argue on the contrary that usual tools and ideas of statistical 
mechanics do apply to such systems, both at equilibrium and out of equilibrium. How- 
ever, reviewing a variety of recent works, we will show that a careful application of 
these tools reveals truly unusual and fascinating behaviors, absent from the world of 
short range interacting systems. 

After a brief introduction on the unusual negative specific heat and other peculiar 
thermodynamical phenomena, we discuss the usual assumptions of equilibrium statis- 
tical mechanics and their interpretation in systems with long range interactions. Based 
on the assumption of equal probability of any configuration with a given energy, we 
then explain why the Boltzmann-Gibbs entropy actually measures the probability to 
observe a given distribution function. This relies on our ability to prove large deviations 
results for such systems. The result of this analysis is that microcanonical and canonical 
ensembles of systems with long range interactions are described by two dual variational 
problems. We explain why such variational problems lead to possible generic ensemble 
inequivalence, and to a richer zoology of phase transitions than in usual systems. A 



natural question then arises: do we know all possible behaviors stemming from long 
range interactions, and, if not, what are the possible phenomenologies? We answer 
this question by discussing a classification of all microcanonical and canonical phase 
transitions, in long range interacting systems, with emphasis on situations of ensemble 
inequivalence n26i\ . Very interestingly many possible phase transitions and situations of 
ensemble inequivalence have not been observed yet. We then describe, for two dimen- 
sional flows, the first observation of appearance of ensemble inequivalence associated 
to bicritical and azeotropy phase transitions. 



Kinetic theories and relaxation toward equilibrium. Because systems with long 
range interactions relax very slowly towards equilibrium, or because they can be forced 
by external field, the study of out of equilibrium situations is physically essential. During 
the past century, there have been many attempts to find a general formalism for out of 
equilibrium statistical mechanics, which would give the equivalent of the Gibbs picture 
for out of equilibrium states. Unfortunately, as recognized by most of the statistical 
mechanics community, until now any such attempt failed. This is mainly due to the fact 
that our knowledge of out of equilibrium situations can not be parameterized by a small 
number of macroscopic quantities, playing the same role as dynamical invariants for the 
equilibrium theory. Then out of equilibrium statistical mechanics must be addressed by 
a case by case careful examination of dynamics, using some appropriate probabilistic 
description. 

For relaxation to equilibrium of Hamiltonian systems with long range interactions, 
standard tools have been developed, mainly kinetic theory. In the introductory para- 
graph, we briefly explain the basic ideas of kinetic theories. In the following, we first 
stress the role of a Vlasov description for small time, and then the role of Lenard-Balescu 
equation (also called coUisional Boltzmann equation in the context of self-gravitating 
systems) for larger time. We also briefly review the recent application of Lynden-Bell 
equilibrium statistical mechanics for the Vlasov equation to simple one dimensional 
models. We also discuss new recent results for the kinetic theory of such systems. 
The first is the generic existence, for the one particle stochastic process, of anomalous 
diffusion and of long relaxation times. We guess that the implications of such a result 
for the validity of the kinetic approach has not been well appreciated up to now. The 
second class of results deals with the time of validity for the Vlasov approximation 
and with the typical time needed to observe relaxation towards equilibrium. One of the 
most striking result is that the Lenard Balescu operator vanishes identically for one 
dimensional systems. This explains the existence of anomalous scaling laws for the 
relaxation towards equilibrium in models like the HMF model. 



Non equilibrium stationary states (NESS). Another class of out of equilibrium prob- 
lems is the study of systems with long range interactions subjected to small non- 
Hamiltonian forces and to weak dissipation. Such a framework is actually the most rele- 
vant one for many physical applications. We will emphasize its interest for geophysical 
flows, like for instance simplified models of ocean currents. 




FIGURE 1. Caloric curve (temperature T = j3 * as a function of the energy U) for the SGR model, 
for £ = 10-2. The energies corresponding to the dashed vertical lines are referred to, from left to right, 
as Uiow, Utup, Uc, and [//„,;/,. The decreasing temperature between Utop and Uc characterizes a range of 
negative specific heat and thus of ensemble inequivalence. At Uc, there is a microcanonical second order 
phase transition associated with the canonical first order phase transition (please see the text for a detailed 
explanation). Such a behavior is linked to the existence of tricritical points in both statistical ensembles 
(see figure |2]i. A classification of all possible routes to ensemble inequivalence is briefly described in 
section 123] or in ll26ll . 



The average balance between forcing and dissipation usually leads to statistically 
stationary states, the properties of which may be studied experimentally, numerically 
and theoretically. As there is no detailed balance, the system is maintained out of 
equilibrium. The fluxes of the Hamiltonian conserved quantities then become essential 
physical variables. 

We show in this last section, that this leads, in the context of two dimensional flows, 
to very interesting out of equilibrium phase transitions. We believe that the study of the 
statistical mechanics of such non equilibrium stationary states and phase transitions, in 
other systems with long range interactions, is one of the main challenges in this field. 



2. EQUILIBRIUM STATISTICAL MECHANICS OF SYSTEMS 
WITH LONG RANGE INTERACTIONS 

2.1. Peculiarities of thermodynamics of systems with long range 

interactions 



For systems with long range interactions, the most intriguing thermodynamical prop- 
erty is the generic occurrence of statistical ensemble inequivalence and negative specific 
heat. Such possibilities have first been recognized and studied in the context of self 
gravitating systems IE livi Ull]. Afterwards, ensemble inequivalence and negative spe- 
cific heat have been observed or predicted in a number of different physical systems: 
two dimensional turbulence EE 12^, plasma physics [|28l. [l9l]. spin systems or toy 
,_or_self_gravitating systems in situations different from the simple initial 
]. A detailed description of each of these cases is provided in 



models [22J25I , 

case fllMi&llo, 
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FIGURE 2. Phase diagrams for the SGR model. Left: for each value of e, the critical energies where 
a microcanonical phase transition occurs are represented by black points or circles. When the softening 
parameter is varied, at the tri critical point e ~ £j, first order phase transitions (black points) change to 
second order phase transitions (circles). Right: the bottom dashed line and black line represent respectively 
Uion- (see fig.[Tl) and C//,,^,/,. The difference C//„i,/, — is the energy jump associated to a first order phase 
transition in the canonical ensemble. This jump decreases to zero at canonical tricritical point Ej where the 
canonical first order phase transition changes to a canonical second order phase transition. The tricritical 
points have different energies and e values in both ensembles. 



To motivate the following development, let us give an example where ensemble in- 
equivalence appears in an unusual way. We discuss the equilibrium properties of the 
Self Gravitating Ring (SGR) model, a toy model for self gravitating systems. Whereas 
we mainly present here its equilibrium properties, we stress that this system is very in- 
teresting also from a dynamical point of view, as it shows a number of out of equilibrium 
quasi-stationary states i3ll [32ll. 

The Hamiltonian of the SGR model is: H = ^Zf-^ pJ - ^if ,-i -4 , 1 . • 

Particles are constrained on a ring (0 < 0, < In). The angles 0/ are conjugate to the 
momenta pi. £ is a small scale softening of the gravitational interaction. We study the 
phase transitions of this system and how they evolve when £ is varied. Please see pOl] 
for a detailed discussion. 

Figure [U shows the caloric curve T{U) for £ = 10^^, where T is the temperature 
and U is the energy. For Utop < U < Uc the temperature is decreasing. The specific 
heat C = dU/dT is thus negative in this area, showing that statistical ensembles are 
not equivalent (in the canonical ensemble, the specific heat is always positive). The 
horizontal dashed line is the Maxwell construction, which links microcanonical and 
canonical ensembles. From it one sees that in the canonical ensemble, when /3 is varied, 
there is a first order phase transition characterized by an energy jump between the values 
Uiow and Uhigh. This is a common feature in case of ensemble inequivalence. 

What is less common is the concomitant existence of a second order phase transition 
in the microcanonical ensemble, at the energy Uc- At this point, the temperature T is 
continuous, whereas its derivative is discontinuous as it is clear from the curve. In figure 
[2l we show that this type of ensemble inequivalence, with the coexistence of a first 
order canonical first order phase transition and of a microcanonical second order phase 
transition, is linked to the existence of a tricritical point in both ensembles. 



The SGR model displays one possible route to ensemble inequivalence, out of sev- 
eral others. However, there are some constraints on the possible phenomenologies. For 
instance, at a second order phase transition, the negative specific heat jump must be pos- 
itive. By contrast, the temperature jumps at a discontinuity associated with a first order 
microcanonical phase transition must be negative (this means that when energy is in- 
creased the system has a negative temperature jump). Summarizing all these constraints 
yields a classification of all possible ensemble inequivalences and their links with 
phase transitions. To prepare the discussion of this classification, we recall now the main 
hypothesis and definitions of statistical equilibrium, in the context of systems with long 
range interactions. 



2.7.7. Additivity extensivity and thermodynamic limit 

When studying statistical mechanics of non additive systems, the first problem one 
has to deal with is the inadequacy of the thermodynamic limit {N ^ °°,V ^ °° -V h 
the volume-, with N /V kept constant). Indeed, what is physically important in order to 
understand the behavior of large systems, is not really to study the large A'' limit, but 
rather to obtain properties that do not depend much on for large A'^ (the equivalent 
of intensive variables). For short range interacting systems, this is achieved through the 
thermodynamic limit; for non additive systems, the scaling limit to be considered is 
different and depends on the problem. Let us consider N particles which dynamics is 
described by the Hamiltonian 

1 ^ iV 

^^ = 2Ei^' + 2 ^ ^(^-^;), (1) 

where c is a coupling constant. The thermodynamic limit in this case amounts to send 
A'^ and the volume to infinity, keeping density and c constant. If V{x) decreases fast 
enough so that interactions for a particle come mainly from the first neighbors, then 
increasing A'^ at constant density has almost no effect on the bulk, and physical properties 
are almost independent of N: the thermodynamic limit is appropriate. This is wrong of 
course if the potential for a particle is dominated by the influence of far away particles. 
The appropriate scaling in this case may be as follows: fixed volume, c oc l/N^, and 
M oo (others equivalent combinations are possible, as the one given below for self 
gravitating particles). 

The best known example of such a special scaling concerns self gravitating stars, for 
which the ration M/R is usually kept constant, where M is the total mass and R is the 
system's radius (thermodynamic limit would be M/R^ constant). Another toy example 
is given and studied for instance in iIbbIi . This type of scaling is also the relevant one for 
point vortices in two dimensional and geophysical turbulence, where the total volume 
and total vorticity have to be kept fixed, but divided in smaller and smaller units. Let 
us note for completeness that in some cases, the thermodynamic limit is appropriate in 



presence of long range interactions, for instance when some screening is involved Il34l1 : 
we shall exclude these cases in the following. 



According to the above discussion, let us rewrite Eq. \T\ using the convenient scaling 

c = ±l/N: 



J 1 ^ 

H^ = -'£pj± V{xt-xj) (2) 

/=1 i,j=l 

This classical scaling of the coupling parameter is called the Kac's prescription (Issl]) 
or sometimes the mean field scaling (see for instance [|36l]). Within this scaling, taking 
the limit N ^ o° with all other parameters fixed (fixed volume for instance), the sum 
over i and j is clearly of order N^, and the energy per particle H^/N is intensive. This 
scaling is also the relevant one in order to obtain the coUisionless Boltzmann equation, 
for the dynamics, in the large limit. We will use Eq. |2]in the following. 



2.7.2. The microcanonical and canonical ensembles 

We suppose that the energy E of our system is known, and consider the microcanon- 
ical ensemble. In this statistical ensemble all phase space configurations with energy E 
have the same probability; the associated microcanonical measure is then 

1 ^ 
LIn[E) 

where Q.^ (E) is the volume of the energy shell in the phase space Q.^ (E) = 
jY[dxidpi5 {Hn {{xi, Pi}) — E). We consider here the energy as the only parameter, 
however generalization of the following discussion to other quantities conserved by the 
dynamics is straightforward. 

The only hypothesis of equilibrium statistical mechanics is that averages with respect 
to will correctly describe the macroscopic behavior of our system. This hypothesis 
is usually verified after a sufficiently long time, when the systems has "relaxed" to 
equilibrium. 

The Boltzmann entropy per particle is defined as 

SN{E) = ^\ogaN{E) 

In the following, we will justify that in the long range thermodynamic limit, the entropy 
per particle S;^ (E) has a limit: 

S^iE) S{E) 

The canonical ensemble is defined similarly, using the canonical measure 

1 ^ 

llc,N = rj WdxidpiQ^^ [-pHf^ iUi^Pi})] , 
^n[E) f-J\ 

with the associated partition function Zp/{P) = jYldridpiCxp [—^H^ {{xi-,Pi})] ^^d 
free energies Fn{[5) = ilogZ^v (/3) and Fa, (j8) F (/3) 



2.2. Large deviation results 



2.2.7. Justification of the Boltzmann-Gibbs entropy 

Let us consider the particle distribution on /i— space / {x,p) (/ {x,p) dxdp is the prob- 
ability to observe a particle with position x and momentum p). f defines a macrostate 
as many microscopic states correspond to a given /. As explained in the previous para- 
graph, the hypothesis of usual statistical mechanics is that all microscopic states with a 
given energy E are equiprobable. Given this uniformity in phase space, we address the 
question: what is the number of microscopic states having the distribution /? 

It is a classical combinatorial result to show that the logarithm of number of micro- 
scopic states corresponding to a distribution / is given by 

s[f] = - J dxdp f log f 

where s is sometimes called the Boltzmann-Gibbs entropy. It is the Boltzmann entropy 
associated to the macrostate /, in the sense that it counts the number of microstates 
corresponding to /. We stress that no other functional has this probabilistic meaning, 
and that this property is independent of the Hamiltonian. 

Thanks to the long range nature of the interaction, for most configurations, the energy 
per particle can be expressed in term of the distribution function /, using 

h{f) ~ ^f{x,p)dxdp + dxidpidx2dp2f{xi^pi)f{x2,P2)V{xi-X2). (3) 

A'^oo J Z J 

This mean field approximation for the energy allows to conclude that the equilibrium 
entropy is given by 

SN{E)=\og{aN{E)) ~ NS{E) with S{E)=sup{s{f) \ h{f) = E} (4) 

In the limit of a large number of particles, the mean field approximation Eq. (|3]) and its 
consequence the variational problem © have been justified rigorously for many systems 
with long range interactions. The first result assumes a smooth potential V and has 
been proved by B6i1 . see also the works by Hertel and Thirring on the self gravitating 
fermions [|27|l. 



2.2.2. Large deviations 



We explained why the Boltzmann-Gibbs entropy is the correct one to describe the 
probability of a given /. Large deviations provide a useful tool to obtain similar results 
in a wider context. We refer to the very interesting contributions of Ellis and coworkers 
( [|37L [29I IbsI [39I S^). We also refer to [i41J for a simple detailed explanation of many 
large deviations results in the context of long range interacting systems. 

In a first step one describes the system at hand by a macroscopic variable; this 
may be a coarse-grained density profile /, a density of charges in plasma physics, a 



magnetization profile for a magnetic model. In the following, we will generically call 
this macroscopic variable m; it may be a scalar, a finite or infinite dimensional variable. 

One then associates a probability to each macrostate m. Large deviation theory 
comes into play to estimate ^{m), the number of microstates corresponding to the 
macrostate m: 

log(tlyv(m)) ~ Nsim) . 

This defines the entropy s{m). 

In a second step, one has to express the constraints (energy or other dynamical 
invariants) as functions of the macroscopic variable m. In general, it is not possible to 
express exactly H; however, for long range interacting systems, one can define a suitable 
approximating mean field functional h{m), as in Eq. ([3]). 

Having now at hand the entropy and energy functionals, one can compute the micro- 
canonical density of states Q.{E) (^J^): the microcanonical solution is simply given by 
the variational problem 

log(f2yv(£)) ~ NS{E) with S{E) = sup{^(m) | h{m) = E} (5) 

In the canonical ensemble, similar considerations lead to the conclusion that the free 
energy and the canonical equilibrium are given by the variational problem 

log(Zyv(/3)) ~ iVF(/3) with F(/3) =inf{-5(m)+/3/z(m)} (6) 

We insist that this reduction of the microcanonical and canonical calculations to the 
variational problems ([5]) and Q is in many cases rigorously justified. 

2.3. Ensemble equivalence and simplification of variational problems 

As discussed in the previous section, the microcanonical and canonical equilibrium 
states are, most of the times, given by ([5]) and Q respectively. These two variational 
problems are dual ones: the canonical one is obtained from the microcanonical one 
by relaxing a constraint. In the following section, we discuss the mathematical links 
between two such dual variational problems. We then apply this to characterize ensemble 
equivalence, and we use it to prove relations between classes of variational problems. 

2.3.1. Relations between constrained and relaxed variational problems 

It is possible to state some general results about the variational problems ^ and Q, 
independently of the precise form of the functions s and h: 

1. a minimizer nic of Q is a minimizer of dS]), with constraint E = h{mc). 

2. a minimizer of ([5]) is a critical point of ^ for some /3, but it is not always a 
minimizer: it is a minimizer of Q if and only if S coincides with its concave hull at 
point E = h{m^). Otherwise, it may be a local minimum, or a saddle point of Q. 



Such results are extremely classical. More detailed results in this context may be 
found in [[btIi . We also refer to [42] for a concise discussion and proof. The previous 
points immediately translate into the language of statistical mechanics, and provide a 
full characterization of ensemble inequivalence: 

• A canonical equilibrium is always a microcanonical equilibrium for some energy E. 

• A microcanonical equilibrium at energy £ is a canonical equilibrium for some 
temperature 1//3 if and only if S coincides with its concave hull at energy E. 
Whenever S coincides with its concave hull, we will say that the ensembles are 
equivalent; otherwise we will say they are not equivalent. 



23.2. Simpler variational problem for statistical equilibria 



In the previous paragraph, we have explained relations between solutions of a con- 
strained variational problem and of the associated relaxed one. Using similar results 
and further theoretical considerations, it is possible to obtain much simpler variational 
problems than the natural microcanonical ones, for the equilibria of Euler and Vlasov 
equations [|42|] . We think that these new results provide essential simplifications that 
will be useful in many studies, we thus describe them in this section. However, from a 
physical point of view, these mathematical results may be viewed as technical, and we 
advise the non expert reader to skip this section at first reading. 



When studying statistical equilibria of systems with long range interactions, one has 
to deal with variational problems with one or several constraints. In the case of the 
statistical mechanics of the Euler (resp. the Vlasov equation), there is actually an infinite 
number of Casimir's functional conservation laws, encoded in the initial distribution d 
of the vorticity field (resp. the particle distribution function). This is a huge practical 
limitation. When faced with real phenomena, physicists can then either give physical 
arguments for a given type of distribution d (modeler approach) or ask whether there 
exists some distribution d with equilibria close to the observed field (inverse problem 
approach). However, in any case the complexity remains: the class of equilibria is huge. 

In the following of this paragraph, we describe recent mathematical results which 
allow to relate the microcanonical equilibria to much simpler variational problems. 
From a physical point of view, this simplification is extremely interesting. We describe 
these results in the context of the equilibrium theory for the Euler equation (Robert 
Sommeria Miller theory ifisi [Till or RSM theory), but the following results may be 
easily generalized to other cases like the statistical mechanics of the Vlasov equation. 
We refer to [;42] for a more detailed discussion. 

>From a mathematical point of view, one has to solve a microcanonical variational 
problem (MVP): maximizing a mixing entropy y[p] = —J^^d^xJ Japlogp, with 
constraints on energy E and vorticity distribution d 



S{Eo,d) = sup {y[p] \E[W]=Eo ,D[p]=d} (MVP). 

{p\N[p]=l} 



p (x, o) depends on space x and vorticity a variables. 

During recent years, authors have proposed ahernative approaches, which led to prac- 
tical and/or mathematical simplifications in the study of such equilibria. As a first exam- 



ple, Ellis, Haven and Turkington [|29|] proposed to treat the vorticity distribution canon- 
ically (in a canonical statistical ensemble). From a physical point of view, a canonical 
ensemble for the vorticity distribution would mean that the system is in equilibrium 
with a bath providing a prior distribution of vorticity. As such a bath does not exist, the 
physically relevant ensemble remains the one based on the dynamics: the microcanon- 
ical one. However, the Ellis -Haven-Turkington approach is extremely interesting as it 
provides a drastic mathematical and practical simplification to the problem of comput- 
ing equilibrium states. A second example, largely popularized by Chavanis [isilU], is 
the maximization of generalized entropies. Both the prior distribution approach of El- 
lis, Haven and Turkington or its generalized thermodynamics interpretation by Chavanis 
lead to a second variational problem: the maximization of Casimir's functional, with 
energy constraint (CVP) 

C(£o,*) =inf I'^Jw] = j^s{(0)d^x \E[(d\=Eq^ (CVP) 

where are Casimir's functionals, and s a convex function (Energy-Casimir functionals 



are used in classical works on nonlinear stability of Euler stationary flows [|45L |46|] . 



and have been used to show the nonlinear stability of some of RSM equilibrium states 



mm). 

Another class of variational problems (SFVP), that involve the stream function only 
(and not the vorticity), has been considered in relation with the RSM theory 

D{G)=mi\l (fx --|Vv/|2 + G(l/A) I (SFVP) 
w y.h [2 J J 

Such (SFVP) functionals have been used to prove the existence of solutions to the equa- 
tion describing critical points of (MVP) ||47|]. Interestingly, for the Quasi-geostrophic 



model, in the limit of small Rossby deformation radius, such a SFVP functional is sim- 
ilar to the Van-Der- Walls Cahn Hilliard model which describes phase coexistence in 



usual thermodynamics 11171. 14811. This physical analogy has been used to make precise 
predictions in order to model Jovian vortices iflTl. l49ll. (SFVP) functionals are much 
more regular than (CVP) functionals and thus also very interesting for mathematical 
purposes. 

When we prescribe appropriate relations between the distribution function d, the 
functions s and G, the three previous variational problems have the same critical points. 
This has been one of the motivations for their use in previous works. However, a 
clear description of the relations between the stability of these critical points is still 
missing (Is a (CVP) minimizer an RSM equilibria? Or does an RSM equilibria minimize 
(CVP)?). This has led to fuzzy discussions in recent papers. Providing an answer is 
a very important theoretical issue because, as explained previously, it leads to deep 
mathematical simplifications and will provide useful physical analogies. 



In Il42ll we establish the relation between these three variational problems. The result is 



that any minimizer (global or local) of (SFVP) minimizes (CVP) and that any minimizer 



of (CVP) is an RSM equilibria. The opposite statements are wrong in general. For 
instance (CVP) minimizers may not minimize (SFVP), but may be instead only saddles. 
Similarly, RSM equilibria may not minimize (CVP) but be only saddles, even if no 
explicit example has yet been exhibited. 

These results have several interesting consequences : 

1 . As the ensemble of (CVP) minimizers is a sub-ensemble of the ensemble of RSM 
equilibria, one can not claim that (CVP) are more relevant for applications than 
RSM equilibria. 

2. The link between (CVP) and RSM equilibria provides a further justification for 
studying (CVP). 

3. Based on statistical mechanics arguments, when looking at the Euler evolution at 
a coarse-grained level, it may be natural to expect the RSM entropy to increase. 
There is however no reason to expect such a property to be true for the Casimir's 
functional. As explained above, it may also happen that entropy extrema be (CVP) 
saddles. 



2.4. Classification of phase transitions 

Beyond the full characterization of ensemble inequivalence we have described above, 
there are many other qualitative features of the thermodynamics that depend only on the 
structure of the variational problems ^ and Q. Indeed, although the precise form of 
the solution obviously depends on the problem at hand through the functions s{m) and 
h{m), it is possible to classify all the different phenomenologies that one may find in 
the study of a particular long range interacting system. The questions in that respect are, 
increasing complexity at each step: 

• what are the different possible types of generic points on an entropy curve S{E) 
(these correspond to different phases)? 

• what are the possible singular points of a generic S{E) curve (these correspond to 
phase transitions)? 

• what are the possible singular points on the S{E) curve, when an external parameter 
is varied in addition to the energy (that is how phase transitions evolve when a 
parameter is varied)? 



We address these different levels in the following paragraphs, using results from pTLllq] . 
These results are obtained by adapting to the dual variational problems ([5]) and ^ 
ideas that lead to the Landau classification of phase transitions. In the long range case 
however, there is no approximation involved, so the classification does not suffer from 
the problems of standard Landau theory (wrong critical exponents for instance). 




cqui iilequil 

inequi2 

FIGURE 3. The three types of generic points; on the left: entropy S{E) curve, on the right: caloric j3 (E) 
curve. Thick, thin and dashed lines correspond respectively to the three types of points. The dotted lines 
shows the Maxwell construction giving the canonical solution in the inequivalence range. 
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FIGURE 4. The three types of phase transitions (Codimension singularities), for a system with no 
symmetry. A: canonical 1 st order transition. B: canonical destabilization (a local minimum of ^ becomes 
a saddle point). C: microcanonical 1st order, temperature jump. 

2.4.1. Generic points of an entropy curve 

There are three types of generic points on the entropy curve, see Fig.[3l 

• Concave points (that is C,, > 0) where canonical and microcanonical ensembles are 
equivalent. 

• Concave points where ensembles are not equivalent. 

• Convex points (Cy < 0), where ensembles are always inequivalent. 

2.4.2. Singular points of a generic entropy curve: phase transitions 

Generic points as described above define segments of entropy curves, separated by 
singular points, that can be of several types. These points for systems without symmetry 
are classified in Fig. IH 
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FIGURE 5. Classification of singularities, when one external parameter is varied (codimension 1 
singularities). The first two columns give the singularity origin and name; the third one gives its status in 
the canonical ensemble: invisible means it has no consequence on the canonical solution; the six following 
plots are entropic and caloric curves showing the crossing of the singularity when the external parameter is 
varied. One recovers the usual phase transitions (triple point, azeotropy, critical point) in both ensembles. 
What is new is the list of the possible links between the behaviors in each ensemble, and the associated 
appearance of ensemble inequivalence. Please see ll26ll for a more detailed explanation, and singularities 
associated with symmetry breaking. 



2.43. Singular points on a singular entropy curve 

When an external parameter is varied, the entropy curve is modified. Some special 
values of the parameter correspond to qualitative changes for the phase transitions. All 
these possible qualitative changes are classified in ll26ll : Fig. [5] summarizes the results. 

All types of phase transitions and ensemble inequivalences found in the literature 
so far are reproduced in the classification. In addition, the classification predicts the 
possibility of new phenomenologies, and new routes to ensemble inequivalence, that 
have not so far been observed in any specific model. 



2.5. Examples of new phase transitions 

In the previous section we explained the existence of many yet unobserved appearance 
of ensemble inequivalence, associated to phase transitions, as described in [|26l] . More 
recently the new finding of two such examples have been reported [[soj]: bicritical 
points (a bifurcation from a first order phase transition towards two second order phase 
transitions) and second order azeotropy (the simultaneous appearance of two second 
order phase transitions at a bifurcation). We present here these new phase transitions; 
they are associated with equilibrium states of the two dimensional Euler equation, when 
there is a linear relation between vorticity and stream function. This was first studied 



by Chavanis and Sommeria 115 111 in the context of the Robert-Sommeria-Miller (RSM) 



statistical mechanics of 2D flows ifisl [Till . They found a criterion for the existence 
of a transition from a monopole to a dipole above a critical energy, for all (closed) 
domain geometry. In this section, we present an alternative method providing the same 
criterion, which generalizes to a large class of models, and thus shows the universality 
of the phenomenon. More interestingly this new method clarifies the nature of the phase 
transitions involved in this problem and makes the link with the existence of an ensemble 
inequivalence region. Those results are presented in a more general context in llsoll . 
where we note the interest of these phase transitions for very simple ocean models. 

Euler equation and associated variational problem. Let us consider the 2D Euler 
equation in a closed domain Si. It can be written as a transport equation for the vorticity 
(0 — AyA: (9f £0 + u.Vft) = 0. The velocity field u is related to O) via the stream function i//: 
u = X Vi//, with i// = on dS. We introduce the projections ft), of the vorticity ft) on a 
complete orthonormal basis of eigenfunctions of the Laplacian: Ae, = A/e,-, where 

the Xi (all negative) are in decreasing order. The stationary states of this equation are 
prescribed by a functional relation ft) = /(i//) . In the following we consider the solutions 
of the variational problem: 

5(£,r) =max{^[ft)] I ^[ft)] = £ &<^[ft)] =r} (7) 
ft) 

The variational problem ([7]) is similar to the generic problem ([5]) studied above, with two 
constraints instead of one. 

• y is the entropy of the vorticity field ft); we restrict ourselves to a quadratic 
functional: [ft)] = -5 («^)^ = -^E w^- 

• (g is the total energy: (g[(o\ = i (^(Vi//)^^ ^ = -jY.i^iCof 

• ^ is the circulation: ^[co] = {0))q, = Y,i i^i) where (e,) = j^ei{x,y)dxdy 

To compute critical points of the variational problem ([7]), we introduce two Lagrange 
parameters /3 and 7, associated respectively with the energy and the circulation conserva- 
tion. Those critical points are stationary solutions for the initial transport equation with 
f{^f) = fi^if — y. The solutions of the variational problem will thus provide the equilib- 
rium states of the Euler equation that present a linear relationship between vorticity and 
stream function, for a given energy and circulation. 

The aim of the following paragraphs is to determine which ones among the critical 
points are solutions of (|7]). It will then be possible to draw a phase diagram in the plane 
{T^E) for those equilibrium states. 

Dual quadratic variational problems. The problem (|7]), with two constraints, will be 
referred to as the microcanonical problem. As already explained earlier, it is sufficient 
to study the easier unconstrained ensembles, unless there is inequivalence of ensembles. 
The strategy is then the following. Start with the easiest problem: 

7(/3,7) = min^{-^[ft)] +/3 <f[ft)]+y<^[ft)]} (grand canonical). Check if aU 
possible values of £ or F correspond to a grand canonical solution; if yes 
the problem is solved, otherwise, we turn to the more constrained problem: 
F(/3,r) = min^ {-^M + /3 <g[q] \ ^[^] = T} (canonical). 



In principle we could eventually have to solve the microcanonical problem. However, 
in this case, we will see that the microcanonical ensemble is equivalent to the canonical 
one: the whole range of E and Y will be covered by canonical solutions. 

We notice first that S are quadratic functionals and that ^ is a linear functional. 

We will thus have to look for the minimum of a quadratic functional with a linear part. 
Let us call Q the purely quadratic part and L the linear part of this functional. Then we 
have three cases 

1. The smallest eigenvalue of Q is strictly positive. The minimum exists and is achieved by a 
unique minimizer. 

2. At least one eigenvalue of Q is strictly negative. There is no minimum. 

3. The smallest eigenvalue of Q is zero (with eigenfunction eo). If Let^ = (case 3a), the 
minimum exists and each state of the neutral direction {aeo} is a minimizer. If Leo ^ 0, 
(case 3b) then no minimum exists. 

The grand canonical ensemble. In that case the quadratic operator Q associated to 
^ = —y + P S" + is diagonal in the Laplacian eigenvector basis. The variational 
problem admits a unique solution if and only if /3 > Ai (case 1. above). If /3 = Ai(case 
3. above), a neutral direction exists if and only if 7=0. By computing the energy 
and circulation of all those states, we prove that there is a unique solution at each 
point in the diagram {E.T), below a parabola ^ (see figures [6] and |7]-a). Because 
the values of energies above the parabola ^ are not reached, we conclude that there 
is ensemble inequivalence for parameters in this region. We then turn to the more 
constrained canonical problem to find solutions in this area. 

The canonical ensemble. The circulation is now fixed. We first transform this prob- 
lem into an unconstrained variational problem. By using the circulation constraint, we 
express one coordinate in term of the others: 0)\ = (F — X^,-ft),(e,)) / (ei). This expres- 
sion is then injected in the functional ^ = —y + /3<f . The problem is now to find a 
minimizer {tt)/},>2 of this functional, without constraints. This case requires more com- 
putations that the previous one since the operator Q associated to the quadratic part of 
^ is no more diagonal in the basis {e,}. 

We first notice that if the domain geometry admits one or more symmetries, it gener- 
ically exists eigenfunctions having the property (e,) = 0. In the subspace spanned by 
all those eigenfunctions, Q is diagonal, and its smallest eigenvalue is positive as long 
as /3 > /3p where /3{^ is the greatest A, on this subspace. Then we look for the value of 
/3 such that the smallest eigenvalue of Q is zero in the subspace spanned by eigenfunc- 
tions with (e,) 7^ 0. Let us call /3* this value, and CO* the corresponding eigenfunction: 
Q[g)*] = 0. We find after some manipulation that /3* is the greatest zero of the func- 
tion f{x) = 1 — xY,i>i{ei)^ / — ^i) ■ We conclude that there is a single solution to the 
variational problem if and only if /3 > max (/3{^, /3*) . When /3 = max (/3p /3*) , we dis- 
tinguish two cases according to the sign of jSj* — j8* to discuss the existence of a neutral 
direction: 

• i) < j8* we then consider j8 = jS*. There is a solution (case 3a) if F = and no solution 
(case 3b) for F^O. 

• ii) > j3* we then consider j3 = Py. There is a solution (case 3a) for all values of P. 
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FIGURE 6. Case i) Phase diagram for a domain having no symmetry, or with a symmetry axis in the 
case p* < jSp a) Internal parameters are circulation F and energy E. Straight line; first order transition. 
Parabola above this line there is no grand canonical solutions; only canonical and microcanonical 
solutions exist. It corresponds to an ensemble inequivalence region, b) examples of flow, c) Equilibrium 
entropy 5 as a function of E and F. We see clearly a region where the concave envelop of 5'(£',r) is 
not equal to 5'(£',r), which includes the first order transition line; d) )3 = dS/dE at fixed F: there is 
no singularity, e) 7= dS/dY at fixed E: there is a discontinuity of 7 at F = (first order transition). 
The Maxwell construction shows that the ensemble inequivalence region is associated with a first order 
transition in the ensemble where the value of the energy is fixed. 



We thus obtain a criterion, namely the sign of /3j — /3*, that provides two classes of 
phase diagrams. This criterion depends only on the domain geometry. If the domain 
admits a symmetry axis, and is sufficiently stretched in a direction perpendicular to this 
axis, then /3{^ > /3*. For an ellipse, this is always the case. However, there are domains 
with a symmetry axis for which one can find a critical aspect ratio Tc that separates the 
two cases. This is for instance the case of rectangular domains. For T > Tc , case ii) is 
realized (/3{' > /3*) and for t < Xc, case i) is realized (/3{^ < /3*). To conclude, we find 
two classes of phase diagrams: 

• i) Domain without symmetry, or with symmetry with aspect ratio T < Tc (see figure 
[6l). If r 7^ 0, all energy values are reached in the canonical ensemble, and there is a 
single solution for each point (r,^). If F = (straight line on figure [61), two states 
coexist for each energy value E, and /3 = /3*. 

• ii) Domain with symmetry and with an aspect ratio T > Tc (see figure |7]). There is a 
unique canonical solution at each point of the diagram below a parabola I^q in the 
plane {E, F) (see the dashed line figure|7]). Above this parabola, /3 = /3{', and there is 



I 



m- ,r....... d) C) 




FIGURE 7. Case ii) Phase diagram for a domain having a symmetry axis, and sufficiently stretched 
in a direction perpendicular to this axis (case J3* > P^). a),b),c),d),e): see figure Dashed parabola 
J^q: second order transition line. Above this line, dS/dE = fif, and the choice of one state among two 
possibilities breaks the system's symmetry. We see that dfi/dE is discontinuous at point G, as well as 
dy/dr at points B and C, which is the signature of a second order phase transition. 



two canonical solutions at each point of the diagram. They differ only by the sign of 
the contribution of the Laplacian eigenvector associated to the eigenvalue (it is a 
dipole). Below the Parabola, this eigenmode has no contribution to the solution. The 
choice of one solution among the two possibilities breaks the system's symmetry 
when tpQ is crossed. At high energy, the contribution of the dipole is dominant. 

In both cases we find that all circulation and energy values have been reached by 
canonical solutions. We conclude that microcanonical and canonical ensembles are 
equivalent: all microcanonical solutions are also canonical solutions. 

Description of phase transitions. The difference between the two classes of flow 
is the existence of either first or second order phase transitions, corresponding respec- 
tively to case i) and ii). The main observation is that first order and second order mi- 
crocanonical transitions always take place in the ensemble inequivalence area. In that 
respect, those transitions are signatures of ensemble inequivalence for this class of two- 
dimensional flows, which could be observed in laboratory experiments on quasi two 
dimensional flows. 

The transition from systems of type i) to systems of type ii), when the geometry 
is modified, leads to one of the predicted, but yet unobserved phase transition of the 
classification [[26|] . described here in section Classification of Phase Transitions. On 



FIGURE 8. Bicritical point (from Ref. ifsoll ') for a rectangular domain with aspect ratio T. Bifurcation 
from a first order transition line to two second order phase transition lines. Insets are schematic projections 
of the Entropy ^[o] in a plane {(0^,(0*), by taking into account the constraints. M: maximum, m: 
minimum, s; saddle. 



figure m we consider a fixed energy, and present the phase diagram in the (F, t) plane, 
where T is an external parameter characterizing the aspect ratio of the domain. In the 
microcanonical ensemble, there is a bifurcation from a first order transition line to two 
second order transition lines at a critical value T = Tc. Such a bifurcation is referred to 
as a bicritical point. 



2.6. Challenges in the equilibrium statistical mechanics 

It seems fair to say that the equilibrium statistical mechanics of long range interacting 
systems is well understood at a fundamental level, despite the important differences with 
the short range case. However, some challenges and questions remain open; we mention 
here some of them. 

A first issue is the relevance for natural phenomenon or laboratory experiments: is 
it possible to identify situations where equilibrium statistical mechanics satisfactorily 



describes the structures observed (like in Jovian atmosphere [Inn )? 

The main challenge is probably to observe negative specific heat and ensemble in- 
equivalence in experiments or in real physical systems. Up to now, this has not been 
possible. A first possibility would be to observe natural phenomena corresponding to sta- 
tistical microcanonical equilibria, which are characterized by negative specific heat, like 
for instance geophysical flows in microcanonical situations. Another possibility would 
be to achieve this in the lab. One should consider a physical system sufficiently simple in 
order to be able to characterize ensemble inequivalence, or the computation of negative 
specific heat. Moreover, exchanges of energy with the environment have to be negligible 
over a sufficiently long time, in order to make sure that the microcanonical ensemble is 
the relevant one; this imposes severe constraints on any laboratory setup. Several can- 
didates have been considered. One of the most popular one may be to build synthetic 



magnetic systems with long range interactions. Another possibility could be to design 
simple two dimensional flow experiments or two dimensional plasma experiments, in 
order to reproduce the recently predicted ensemble inequivalence, as briefly describe 
inQ. 

There also several theoretical issues : 

• Could we find physical examples of the new phase transitions found in the classifi- 
cation? 

• The structure of the dual variational problems ^ and ^ appears in other physical 
contexts. Could the results described here, like for instance the classification of 
phase transitions, have an interest, when applied to these different situations? For a 



step in this direction, see for instance 115211 . 



• A very interesting and difficult challenge would be to make the classification 112611 
rigorous. This implies to give a precise mathematical definition to the notion of a 
normal form for a mean field variational problem. 

3. KINETIC THEORIES OF SYSTEMS WITH LONG RANGE 

INTERACTIONS 

The previous section provided a brief summary of old and new results for long range 
interacting systems at equilibrium. Unfortunately, it turns out that the coherent structures 
these systems form, and the stationary states they reach are generally out of equilibrium. 
Although knowledge of equilibrium is a useful benchmark, which usually yields a 
qualitative understanding of the physics, some new techniques are needed to really 
understand the phenomena at hand. Clearly, we have to reintroduce the time in our 
framework and study the dynamics of the systems. 

We have seen that for systems with long range interactions, a mean field approach is 
usually exact in the limit of a large number of particles, when one wants to describe the 
equilibrium macrostates. This is valid thanks to an averaging of the potential over many 
particles. In the following we explain that a similar mean field approach is also valid 
for the dynamics: at each time the potential and the force can be expressed with a very 
good approximation from the one particle distribution function, and thus the BBGKY 
hierarchy can be safely truncated. 

This well understood fact is the base of the kinetic theory for the dynamics of 
systems with long range interaction. This led to the classical kinetic theories of self 
gravitating stars, plasmas in the weak coupling limit, or point vortex models in two 
dimensional turbulence. In the limit of a large n umber of particles, such dynamics is 
well approximated by kinetic theories ll53l.l54ll6l.l55l.l56l.l57n: to leading order in 1/ y/N 



the dynamics is of a Vlasov type; after a much longer time, the relaxation towards 
equilibrium is governed by Lenard-Balescu type dynamics (or its approximation by the 
Landau equation). 

In the next subsection, we introduce briefly the Vlasov dynamics, and the issue of 
its time of validity. For a large number of particles, these systems may exhibit quasi- 



stationary states (QSS) ||58L l59 |l (in the plasma or astrophysical context see for instance 



ll60Ll53|]). We give a kinetic interpretation of such states as stable stationary solutions of 



the Vlasov dynamics. 

An interesting question is whether we can predict such Quasi-Stationary States, from 
the initial condition of the Vlasov equation, using statistical mechanics. In the following 
subsection, we present recent studies on the equilibrium statistical mechanics of the 
Vlasov equation (and not of the N particle dynamics) 116 ll . l62l . 16311 . in the spirit of Lynden- 
Bell's work (jij] in the context of self-gravitating stars. 

We then turn to the kinetic theory of these systems beyond the time of validity of 
Vlasov equation, which leads to the Lenard-Balescu equation. This allows to address 
the important question of the time scale for the relaxation to equilibrium: this time scale 
may be of order N/logN (this is a classical result by Chandrasekhar for relaxation to 
equilibrium of self-gravitating stars, or of a plasma), of order A'^ (for a smooth potential), 
or much larger than A'^ (this is related to the recent result that the Lenard-Balescu operator 
identically vanishes for one dimensional systems fs^j). This last result explains the 
striking numerical observation of an A^^^ -^ time scale in the HMF model [i59l1 . 

In the following subsection, we explain how a classical kinetic approach allows to 
describe the stochastic process of a single particle in a bath composed by a large number 
of other particles. This stochastic process is governed by a usual Fokker-Planck equation. 
In classical papers, this bath is at equilibrium. We stress here that this bath can also be a 
bath of particles in an out of equilibrium Quasi Stationary State. We explain recent new 
results 115 511 proving that this Fokker Planck equation has no spectral gap, and lead to 
long time algebraic correlations and anomalous diffusion. This provides a quantitative 
prediction for the algebraic autocorrelation function and anomalous diffusion indices, 
previously observed in some numerical computations f,64i . ,65^ ,661 . These theoretical 
predictions have been numerically checked in [67]. Some more recent related results 
have also been reported in (6§]. We note that an alternative explanation, both for the 
existence of QSS and for anomalous diffusion has been proposed in the context of 



Tsallis non extensive statistical mechanics |i69|, |66i] (see 115911 and 1551, lZOQ for further 
discussions). 

The last subsection is devoted to describe some remaining issues and challenges in 
the context of the classical kinetic theory for systems with long range interactions. 
We also note that we do not describe many other existing dynamical properties which 
are common to systems with long range interactions: vanishing Lyapounov exponents 
[SSLIvIll, breaking of ergodicity [T^, Ivj, [vij, ans so on. All of the common dynamical 
properties of systems with long range interactions are a result of similar collective (self- 
consistent) dynamics [|75|l . 



3.1. Vlasov dynamics and Quasi-Stationary states 

As for the equilibrium statistical mechanics, one needs to choose a scaling to study 
the kinetic theory; again, the scaling described in the equilibrium context, which ensures 
that each particle experiences a force of order 1, is the appropriate one. The goal is now 
to approximate the dynamics of N ordinary differential equations for the discrete parti- 
cles dynamics by a single Partial differential equation for the one-particle distribution 
function. 



For definiteness, we consider the following Hamiltonian system: 



/ ^' = (9^ 

t Pi - -Nl^m^K^i-^j) 

The range of the potential V is supposed to be of the same order of magnitude as the 
total size of the system: this is our definition for "long range interaction' 
Consider the following continuous approximation of the potential: 



^{x,t) = jv{y-x)f{y,p,t)dydp; 



(9) 



and the corresponding equation for the one-particle distribution function / (this is the 
Vlasov equation): 

^+p^_^^=0 (10) 
dt dx dx dp 

Replacing the true discrete potential by 4> neglects correlations between particles and 
finite-A^ effects. However, as each particle interacts at any time with an extensive number 
of other particles, one may hope that this mean field approach correctly reproduces the 
potential experienced by a particle, and becomes exact in the infinite N limit. Under 
some regularity assumptions for the potential V , this is indeed correct, and it has been 



rigorously proved (see [|76|] for a very regular V, [|77|] for a mildly singular potential). 
To be more precise, these theorems state the following: take a discrete A^-particles initial 
condition and an initial continuous one-particle distribution function /(x,p,0) which is 
close, in some sense, to the former; evolve the A^^ particles according to ([8]), and evolve 
/(jc,p,0) according to ^ and (flOl) : then the A/^-particles dynamics and f{x,p,t) will 
remain close for a time at least of the order of logA'^. Several remarks are in order: 

1. This implies that if the f ^ oo limit is taken for a fixed A'^, finite-A'^ effects will come 
into play; the evolution will then depart from the Vlasov equation, and we expect 
the system to eventually reach the statistical equilibrium. However, for any finite 
time T there exists some A'^ such that the system approximately follows the Vlasov 
equation up to time T . 

2. The logA'' is optimal in the sense that there exist initial conditions such that the 
discrete ^ and Vlasov (flOl) dynamics diverge on such a time scale (see [78] for 
further discussion). H 

3. However, this "coincidence time" may in some cases be much longer: for instance, 
discrete initial conditions close to a stable stationary state of the Vlasov equation 
stay so for a time algebraic in A^ (see [|59|] for a numerical observation and [[80|] for 
a mathematical investigation of the phenomenon). 



^ Albeit rather general, equations (|8]i do not include the 2D flows, nor the wave-particles models; most of 
the following discussion does apply to these cases too, with small modifications. 

^ A recent consideration of the thermodynamic stability of a mean field Ising model with stochastic 
dynamics has found the relaxation time to be logarithmic in ll79ll . 



4. The analogous result for 2D flows is the convergence of the dynamics of discrete 
vortices to the corresponding continuous partial differential equation for the vor- 
ticity field (Euler, Quasi-geostrophic...)- For 2D flows however, the fundamental 
equation is the PDE, contrary to the classical particles case. A mathematical proof 
of convergence is given in jSlll . For wave-particles systems, the analogous theorem 
is given in [82]. 

5. The mathematical proofs cited above do not include the gravitational and electro- 
static l/r cases. It seems however reasonable to believe that some convergence 
result towards the Vlasov equation still holds in this case; Vlasov equation is rou- 
tinely used by physicists for these potentials. 

In the light of the previous remarks, and if the number of particles is big enough, the 
following dynamical scenario now seems reasonable: 

• Starting from some initial condition, the A/^-particles system approximately follows 
the Vlasov dynamics, and evolves on a time scale of order 1 . 

• It then approaches a stable stationary state of the Vlasov equation; the Vlasov 
evolution stops. 

• Because of discreteness effects, the system evolves on a time scale of order A^" 
for some a, and slowly approaches the full statistical equilibrium, moving along a 
series of stable stationary states of the Vlasov equation. 

In this scenario, the A/^-particles system gets trapped for long times out of equilibrium, 
close to stable stationary states of the Vlasov equation: these are then called "quasi 
stationary states" in the literature. This is the basis for the "violent relaxation" theory of 
Lynden-Bell fl*]; Refs. [61, 59, 83] give examples of this scenario for wave-particles, the 



HMF and astrophysical models. The next problem is then to study the stable stationary 
states of the Vlasov equation, that is the candidates for the "quasi stationary states". 
Before turning to this in the next paragraph, let us note that there is however no reason for 
this scenario to be the only possibility: for instance, the Vlasov dynamics may converge 



to a periodic solution of the Vlasov equation 118411 . 

The Vlasov equation (as well as the Euler equation and its variants) has many in- 
variants: beside the energy H[f] (and possibly the linear or angular momentum), inher- 
ited from the discrete Hamiltonian equations, the following quantities Cs[f], sometimes 
called Casimirs, are conserved /or any function s: 

Cs[f]= js{f{x,p,t))dxdp. (11) 

Using these invariants, it is possible to construct many stationary states of the Vlasov 
equation. Consider the following variational problem, for a concave function s: 

sup I s{f{x,p)) dx dp I j fdxdp=l , //[/]= e| . (12) 



Any solution of this variational problem yields a stationary solution of the Vlasov 
equation. In addition, the variational structure of the construction is very useful to study 
the stability of such states (see [2§,[85l] for more details). There is no constraint on the 



concave function s, so that we have very many stable stationary states of the Vlasov 
equation. As a consequence, many numerical or experimental results can be fitted with a 
good choice of s; this is also a serious limit of the theory: without a recipe to choose the 
right stationary state, the theory is not predictive. We address this problem in the next 
paragraph. 

We have explained that any Vlasov stable stationary solution is a Quasi Stationary 
State. Then, because inhomogeneous Vlasov stationary states do exist, one should not 
expect Quasi Stationary States to be homogeneous. This is illustrated in the case of 
several generalizations of the HMF model in Ref. 117 811 . 

The issue of the robustness of QSS when the Hamiltonian is perturbed by short range 
interactions [,861] or when the system is coupled to an external bath [|87[| has also been 
addressed, and it was found that while the power law behavior survives, the exponent 
may not be universal. 



3.2. Lynden Bell statistical mechanics 

Under the Vlasov dynamics, the distribution function / is advected, by a field which 
itself depends on /. The conservation of Casimirs amounts to the conservation of the 
area of all level sets 

I[a,b] = ii^^P) su'^h that a < f{x,p,t) < b} ; 

as time goes by the sets are filamented down to a finer and finer scale, and the filaments 
get interwoven. Understanding the long time behavior of this complicated dynamics is 
not an easy task, analytically or numerically. Assuming that / tends to one of the many 
stable stationary states of the Vlasov equation, the Lynden-Bell statistical mechanics 
is a recipe to choose the right one. In a nutshell, at fixed energy (and possibly linear 
or angular momentum), it selects the most mixed state compatible with the Casimirs 
conservation. It is a maximum entropy theory; the Lynden-Bell equilibrium is given by 
the solution of a problem like (fT2)) . the function s being determined by probability theory 
and the initial distribution f{x, p,t = 0) E|. 

The idea goes back to a pioneering work of Lynden-Bell in the context of astrophysics 
in 1967 [1]; the problem was later revisited by Chavanis and collaborators [|88D, in 
connexion with the statistical mechanics of 2D flows. Let us note that the analog of 
Lynden-Bell theory for the Euler and Euler-like equations of 2D flows is the Robert- 
Sommeria-Miller theory [89, Hi [13]: it relies on the very same ideas. 



The Lynden-Bell and Robert-Sommeria-Miller theories have had important successes; 
let us mention here the descriptions of the core of elliptical galaxies, and the giant 
vortices in Jupiter's atmosphere (iT^ . However, this is the exception rather than the 
rule. A lot of works have been devoted to checking these theories in different contexts, 
to which we do not do justice here. To summarize them very briefly, the rule is that 



^ We have to mention that the determination of the Lynden-Bell equilibrium is in general a difficult task; 
the calculations are usually practical only for two- or three-levels initial distributions. 



the phase space mixing induced by the Vlasov equation is not strong enough, so that 
the theoretical predictions are in general at best qualitatively correct (see ll63n for a 
discussion of these issues; see also [i9C^1 ). 



3.3. Order parameter fluctuations and Lenard-Balescu equation 

In this section, we explain briefly how one classically obtains exact expressions for 
the 1 / fluctuations of the order parameter, for a system with long range interactions 
close to a Quasi Stationary State. In order to make this discussion as simple as possible, 
we treat the case of the HMF model, a one dimensional system with a smooth two 



bod y po tential V. We follow pSi] . and refer to [|91l] for a plasma physics treatment, to 
[56, 57, 1^ for the case of point vortices and to Ref. [92] for self-gravitating stars. 

One could use an asymptotic expansion of the BBGKY hierarchy, where 1 / ^/N is the 
small parameter, and obtain the same results. The \/^/N fluctuations would then have 
been obtained by explicitly solving the dynamical equation for the two point correlation 
function, while truncating the BBGKY hierarchy by assuming a Gaussian closure for 
the three point correlation function. Such a procedure is justified in the large limit 



(see Ref. [|54|]). Our presentation rather follows the Klimontovich approach. 

The state of the A^^-particles system can be described by the discrete single particle 
time-dependent density function (r,jc,j!7) = jjY!j=iS ^ {p~Pj{^)) •, where 

d is the Dirac function, {x,p) the Eulerian coordinates of the phase space and (xi.pi) 
the Lagrangian coordinates of the particles. The dynamics is thus described by the 
Klimontovich 's equation jsill . 

^+ ^_^^ = (13) 
dt dx dx dp 

where the potential V that affects all particles is V{t,x) = — jQ^dy J^^dp cos(x — 
y) fdit^y^p)- This description of the Hamiltonian dynamics derived from ^ is exact : as 
the distribution is a sum of Dirac functions it contains the information on the position 
and velocity of all the particles. It is however too precise for usual physical quantities 
of interest but will be a key starting point for the derivation of approximate equations, 
valid in the large A'^ limit and describing average quantities. 

When N is large, it is natural to approximate the discrete density by a continu- 
ous one f{t,x,p). Considering an ensemble of microscopic initial conditions close to 
the same initial macroscopic state, one defines the statistical average (fd) = fo{x,p), 
whereas fluctuations of probabilistic properties are of order 1 / ^/N. We will assume that 
/o is any stable stationary solution of the Vlasov equation. The discrete time-dependent 
density function can thus be rewritten as fd{t,x,p) = foix^p) + 5f{t,x,p) / \/N, where 
the fluctuation 5f is of zero average. We define similarly the averaged potential {V) and 
its corresponding fluctuations 5V{t,x) so that V{t,x) = {V) + 5V{t,x)/ VN. Inserting 
both expressions in Klimontovich 's equation ([T3]) and taking the average, one obtains 

d^^^dfy_d(y)_dJo ^ ]_/d5Vd5f\ ^^^^ 
dt dx dx dp N \ dx dp / 



The Ihs is the Vlasov equation. The exact kinetic equation (fT4l) suggests that the quasi- 
stationary states of sections [BTI and [X2l do not evolve on time scales much smaller than 
A^; this would explain the extremely slow relaxation of the system towards the statistical 
equilibrium. 

Let us now concentrate on stable homogeneous distributions fo{p), which are station- 
ary since (V) = 0. Subtracting Eq. (fT4l) from Eq. (fT3l) and using fd = fo + Sf/ \/N, one 
gets 



d5f ^ ddf _ d5Vdfo_ 1 
dt dx dx dp y/N 



ddVddf / ddVddf 
dx dp \ dx dp 

For times much shorter than we may drop the rhs encompassing quadratic terms 
in the fluctuations. The fluctuating part 5f are then described, by the linearized Vlasov 
equation (this is another result of the Braun and Hepp theorem 93|]). This suggests 



to introduce the spatio-temporal Fourier- Laplace transform of 5f and dV. This leads to 
dV{0},k) = — / dp—— -, (15) 

where 

dfo 

dp j^^^^r^ (16) 

is the dielectric permittivity. The evolution of the potential autocorrelation, can therefore 
be determined. For homogeneous states, by symmetry, {5V{(0\,ki)5V{Gyi,k2)) = 
except if k\ = —^2 = ±1. 



3.3.1. Autocorrelation of the potential 
One gets, after a transitory exponential decay, the general result 

{5V{n,±l)5V{t2,Tl)) = U dcoe-'^^'^^'^^ 7^^- d^) 



This is an exact result, no approximation has been done yet. 



3.3.2. Lenard Balescu equation 

A similar, but longer, calculation allows to compute the rhs. of Eq. (fT4l) . at order \/N. 
This is very interesting as it gives access to the slow evolution of the distribution /o 
due to the "coUisional" effects. This is, for systems with long range interactions, the 
analogue of the Boltzmann equation for dilute system with short range interactions. We 
do not describe the computation in details (see [19 iL 15411 ). as we just want to discuss 



qualitatively the collision operator. This collision operator is called the Lenard Balescu 
operator and it leads to the Lenard Balescu equation. 

For system of particle with long range interactions given by a two body potential 
xi — X2), the Lenard Balescu equation reads : 



dt 



]_d_ 



dkdp'- 



m 



dfo 



k. /o(p)3^(pO-/o(p')3^(p) 



^p 



5(k. (p- 
(18) 

where k is a wave vector, ^{k) is the Fourier transform of the potential V{x), and 
|£(^, k.p')| is the dielectric permittivity. One note that this is a quadratic operator, as 
for the Boltzmann equation. Moreover, this operator involve a resonance condition in 
the Dirac distribution d (k. (p — p')) • 

>From this equation one clearly expects a relaxation towards equilibrium of any 
Quasi-Stationary state with a characteristic time of order A^. We note that for plasma 
or self gravitating systems, due to the small r divergence of the interaction potential, 
the Lenard Balescu operator diverges at small scales. This is regularized by close two 
body encounters, fixing a small scale cutoff. This leads to a logarithmic correction to the 
relaxation time, which is then the Chandrasekhar time proportional to \og{N) /N. 

One clearly sees on equation (fTSi) that the mechanism for evolution of the distribution 
function is related to the resonances of two particles. An essential point is that the 
condition k. (p — p') =0 cannot be fulfilled for one dimensionnal systems. It would 
indeed imply p = p', and because the Lenard Balescu operator is odd in the variable 
p, it will vanish. Another way to obtain the same result, is to directly compute the rhs 
of Eq. (fT4l) . We do not report such long and tedious computations, but it shows that it 
identically vanishes at order 1 /N, for one dimensional systems. 

This proves that Vlasov stable distribution function will not evolve on time scales 
smaller or equal to A'^. This is an important result: generic out of equilibrium distribu- 
tions, for one dimensionnal systems, evolve on time scales much larger than N. This is 
in agreement with the A'^''^ scaling law which was numerically reported itsT 
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FIGURE 9. Panel (a) presents the temporal evolution of the magnetization M(f ), for the HMF model, 
for different particles numbers: = 102(10^), 10^(102), 2.10^(8), 5.10^(8), 10'^(8) and 2.10'*(4) from 
left to right, the number between brackets corresponding to the number of samples. The horizontal line 
represents the equilibrium value of M. Panel (b) shows the logarithmic timescale ^(A^) as a function of A^, 
whereas the dashed hne represents the law lO'''^^ ^ N^'"^ . From Ref. lf59ll . 



3.3.3. The stochastic process of a single particle in a bath 



Let us now consider relaxation properties of a test-particle, indexed by 1, surrounded 
by a background system of (A'' — 1) particles with a homogeneous distribution. The 
fluctuation of the potential is thus 

5V{t,x) = - dy dp cos{x- y) 5 f{t,y,p) y=cos{x-xi). (19) 

Jo J-oo ^JN 



Using the equations of motion of the test particle and omitting the index 1 for the sake 
of simplicity, one obtains p{t) = p{0) — Jlflu{d5V {u,x{u))) / (dx) / y/N . By introducing 
iteratively the expression of x in the rhs and expanding the derivative of the potential, one 
gets the result at order 1/A^. The key point is that this approach does not use the usual 
ballistic approximation. As a consequence, we obtain an exact result at order l/N. This 
is of paramount importance here to treat accurately the collective effects. As the changes 
in the impulsion are small (of order 1 / v^), the description of the impulsion stochastic 
process by a Fokker-Planck equation is valid. This last equation is then characterized by 
the time behavior of the first two moments {{p{t) — p{0))"). Using the generalization of 
formula (flTl) when the effect of the test particle is taken into account, one obtains in the 
large f -limit 

((p(')-pm ™ 

{{p{t)-p{0)f) ~ (21) 
where the diffusion coefficient D{p) can be written as 

?rVo(p ) 

These results are the exact leading order terms in an expansion where 1/A^ is the small 
parameter. 

Using time variable T = t/N as suggested by Eqs. (l20l) and (|2TI) . the Fokker-Planck 
equation describing the time evolution of the distribution of the test particle is 



D(p) =2Re rdte'P' {dV{t,\)dV{0,-\)) = I!Me}^, (22) 
Jo \£(d.I)\ 



dr dp 



(23) 



We stress that this equation depends on the bath distribution /q. It is valid both for equi- 
librium and and out of equilibrium /o, provided that /o is a stable stationary solution 
of the Vlasov equation. In the limit t ^ oo (more precisely 1 << T << A''), the bracket 
vanishes: the PDF f\ of the test particle converges toward the quasi-stationary distribu- 
tion /o of the surrounding bath. This is in complete agreement with the result that /o is 
stationary for time scales of order A'^. 

All the results of this section, except the fact that the Lenard Balescu equation 
vanishes for one dimensional systems, are classical results. In the next section we explain 
recent results related to the very interesting and peculiar properties of the Fokker-Planck 
equation (l23l) . 



3.4. Autocorrelation function with algebraic decay and anomalous 

diffusion 



In this subsection, we present recent results [|55|1 which predicted the existence of 
non exponential relaxation, autocorrelation of the momentum p with algebraic decay at 
large time, and anomalous diffusion of the spatial or angular variable x. They clarify the 
highly debated disagreement between different numerical simulations reporting either 



anomalous [64] or normal [|65|] diffusion, in particular by delimiting the time regime for 
which such anomalous behavior should occur. We briefly recall that when the moment 
of order n of the distribution scales like t"/^ at large time, such a transport is called 



normal. However, anomalous transport 11941 19511. where moments do not scale as in the 



diffusive case, were reported in some stochastic models, in continuous time random 
walks (Levy walks), and for systems with a lack of stationarity of the corresponding 
stochastic process 1196ll . 

These results have been obtained by analyzing theoretically the properties of the 
Fokker-Planck equation (|23|) . >From the physical point of view, as particles with large 
momentum p fly very fast in comparison to the typical time scales of the fluctuations 
of the potential, they experience a very weak diffusion and thus maintain their large 
momentum during a very long time (one sees from equation (l22j) . using |e(p, 1) p 1, 

that the diffusion coefficient decays as fast as the bath distribution /o {p) for large times). 
Because of this very weak diffusion for large p, the distribution of waiting time for 
passing from a large value of p to a typical value of p, is a thick distribution. This 
explains the algebraic asymptotic for the correlation function. From a mathematical 
point of view, these behaviors are linked to the fact that the Fokker-Planck equation (l23l) 
has a continuous spectrum down to its ground state (without gap). This leads to a non 
exponential relaxation of the different quantities and to long-range temporal correlations 



H55L 19711. These results will generalize to the kinetic theory of any system for which the 
slow variable (here the momentum) live in an infinite space. 

By explicitly deriving an asymptotic expansion of the eigenvalues and eigenfunctions 
of the Fokker Planck equation, the exponent for the algebraic tail of the autocorrelation 



function of momenta has been theoretically computed ||55l . 1971] . This mechanism is 
new in the context of kinetic theory. However, we have discovered later that similar 
Fokker-Planck equations, with a rapidly vanishin g d i ffusi on coefficients obtained by 



other physical mechanisms, had been studied |98i 9^ 100 1. A more recent alternative 



upon. 

161, 



approach to the same phenomena has been proposed [681], together with interesting 
discussions of kinetic applications. 

Let us present the results in the context of the HMF model, for which algebraic large 
time behaviors for momentum autocorrelations had been first numerically observed 
in Refs. ||69l [66|] . In its Quasi Stationary States, the theoretical law for the diffusion 



of angles cr^(T) has been also derived in ll55l . I97ll . The predictions for the diffusion 
properties are listed in Table [B 

When the distribution /o(p) is changed within the HMF model, a transition between 
weak anomalous diffusion (normal diffusion with logarithmic corrections) and strong 



TABLE 1. Asymptotic forms of initial distributions foip), and theoret- 
ical predictions of correlation functions Cp{T) and the diffusion (J,'(t) in 
the long-time regime. Asymptotic forms of the distribution and the pre- 
dictions are assumed and predicted in the limits |/?| ^ °o and t ^ oo re- 
spectively, where T = t/N is a rescaled time. The exponent a is given as 
a = (v-3)/(v + 2). SeeRef. iHlsl for details. 
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FIGURE 10. Angles diffusion ( < {x{t) -x(0))^ > as a function of time) in the HMF model, for a 
quasi-stationary state. Points are from a N body numerical simulation, the straight line is the analytic 
prediction by the kinetic theory. For large times < (x(f) —x{Oj)^ > ^ with V 7^ 1. Such an unexpected 

anomalous diffusion is also observed at equilibrium (see ifssll for more details) 



anomalous diffusion is thus predicted. We have numerically confirmed the theoretical 
predictions [67]. For initial distributions with power-law or Gaussian tails, correlation 
function and diffusion are in good agreement with numerical results. Diffusion is indeed 
anomalous super-diffusion in the case of power-law tails, while normal when Gaussian. 
In the latter case, the system is at equilibrium, but the diffusion exponent shows a loga- 
rithmically slow convergence to unity due to a logarithmic correction of the correlation 
function. This long transient time to observe normal diffusion, even for Gaussian dis- 
tribution and at equilibrium, suggests that one should be very careful to decide whether 
diffusion is anomalous or not. 

We note the existence of another interpretation of Quasi Stationary States and anoma- 
lous diffusion : the algebraic behaviors for momentum autocorrelations have been fitted 
using q-exponential functions i69l [66ll. derived from Tsallis' non extensive statistical 
mechanics. Our theoretical and numerical results are in disagreement with this interpre- 
tation (see ll73,|55i[67l] for further discussions). By contrast with the use of non extensive 



statistical mechanics, we think that our theory explains the phenomena of long range 
temporal correlation and of anomalous diffusion from first principles. 



3.5. Challenges in kinetic theories 

Our main message in this section is that a classical kinetic theory approach for these 
long range interacting systems already explains many intriguing behaviors of these 
systems. However, in contrast with the equilibrium theory, many questions remain open; 
we mention here some of them, without any pretention to be exhaustive: 

• Can we find a better recipe than Lynden-Bell's theory to predict the outcome of 
the Vlasov evolution? This seems hopeless in a general setting (see for instance the 
discussion in ll63ll ). 

• Is it possible to explain the 1.7 exponent for the relaxation to equilibrium in the 
HMF model, and does it have some universality? More generally, is it possible to 
extract other general features of the dynamics beyond the Vlasov equation, like the 
anomalous diffusion, or the long relaxation times described above? 

• At the mathematical level: is it possible to improve on [80] concerning the lifetime 
of QSS? Can the convergence theorems to the Vlasov equation be extended to more 
singular potentials? 

• The most important issue concerning kinetic theories is a clear understanding of 
the limits of validity the different equations. Whereas, for smaller times, kinetic 
theory are based on solid theoretical arguments, the understanding of larger time 
behavior of an ensemble of trajectories, initially close to one another, is not yet 
understood. Numerical computations could be very useful in order to understand 
that. Very few direct numerical tests of the kinetic theories have been performed 
up to now. The main reason is probably the difficulty for such tests, because of the 
long time needed for such test. We think it would be highly relevant to consider 
such problem, in models as simple as possible. 

4. OUT OF EQUILIBRIUM 
4.1. Motivations 

We have described the computation of equilibrium states for systems with long range 
interactions in the first section, and addressed the problem of relaxation to equilibrium 
in the second one. These two types of problems concern isolated Hamiltonian systems, 
systems which may be considered so on the relevant time scales, or systems in contact 
with a thermal bath. In many cases of interest, the system experiences random forces 
and dissipation. Very often the mechanism for dissipation and random forces are from a 
different origin, and do not act as a thermal bath. As a consequence, detailed balance is 
no more valid and the system is subject to fluxes of energy or possibly of other conserved 
quantities; the average energy of the system is fixed by the balance between forcing and 
dissipation. The understanding of the properties of the corresponding Non Equilibrium 



Steady States (NESS) is thus of deep importance. We present here first studies of such 
NESS in the context of systems with long range interactions. The most prominent result 
is the finding of out of equilibrium phase transitions. 

These first studies have been done in the context of two dimensional flows. This 
is indeed essential in this case, as in many applications of fluid dynamics, one of the 
most important problem is the prediction of the very high Reynolds' large-scale flows. 
The highly turbulent nature of such flows, for instance ocean circulation or atmosphere 
dynamics, renders a probabilistic description desirable, if not necessary. At equilibrium, 
a statistical mechanics explanation of the self-organization of geophysical flows has been 
proposed by Robert-Sommeria and Miller (RSM). Out of equilibrium, there are several 
practical and fundamental problems to understand: How the invariants are selected by 
the presence of weak forces and dissipation? What are the associated fluctuations? Are 
all forcings compatible with RSM equilibria? 

We will thus study the Navier Stokes equation with weak random stochastic forces 
and dissipation: 

d(0 

— +u.Vco = vAco-aco + fs (24) 
ot 

where CO is the vorticity, is a random force, act) is the Rayleigh dissipation and v is 
the fluid viscosity. 



4.2. Out of equilibrium phase transitions 

In many turbulent geophysical flows, one can see transitions, at random times, be- 
tween two states with different large scale flows. The most famous example are probably 
the time reversal of the earth magnetic field. We may cite also experimental studies of 



such phenomena, for two dimensional magneti c flow s [|101|] . rotating tank experiment in 



relation with weather regimes in meteorology lll02ll . or magnetic field reversal in MHD 



[|103|] . In all these examples, this generic phenomenon takes place in systems with a 
large number of degrees of freedom. The case of simple turbulent flows may be stud- 
ied in much details theoretically and numerically; we focus here on the case of the two 
dimensional Navier-Stokes equation with a random force. 

Figure \TT\ shows the relaxation of the 2D Navier-Stokes equation to a statistically 
stationary state. It illustrates that, depending on the aspect ratio of the domain, two types 
of large scale flows are possibly observed, either dipoles or unidirectional (zonal) flows. 
We note that these two topologies are also predicted by the equilibrium statistical theory. 

As shown on figure [12l for some values of the control parameter (the aspect ratio 
of the domain), we observe the coexistence of these two flow topologies: the system 
switches back and forth, at random times, between dipole and unidirectional flows. This 
phenomenology is similar to what happens when noise is added to a bistable system. A 
crucial difference here, is that the deterministic dynamics does not have two different 
attracting states (there is no double well potential in this case). 

These few figures show that NESS, for systems with long range interactions, may 
exhibit very interesting phenomena. We hope that this will open a large number of new 
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FIGURE 11. Right panel: relaxation towards a statistically stationary state of the two dimensional 
Navier Stokes equation. After a transient state with a dipole vorticity field, the flow switches to a zonal 
(unidirectional) organization of the vorticity field. Left panel: for a different value of the control parameter 
(here the aspect ratio of the domain), we observe a dipole organization in a statistically stationary situation. 
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FIGURE 12. Right panel: \zi |, the modulus of the first Fourier mode of the vorticity field in the direction 
parallel to the unidirectional flow (for unidirectional flows, \zi \ is close to 0, while for dipoles it oscillates 
around 0.5. The flow thus switches randomly from dipole to unidirectional flows. Left panel: PDF of the 
complex variable zi- 



fundamental works on the subject. Moreover, in a forthcoming paper, we will discuss 
the application of these results to geophysical flows. 

5. CONCLUSION 

We have briefly reviewed in this contribution old and new results on the old, but active, 
subject of systems with long range interactions. We clearly acknowledge that this review 
is far from exhaustive and represents our personal interests. 

In conclusion, it seems to us that equilibrium statistical mechanics of these non addi- 
tive systems is very well understood: a careful application of standard tools allows one 
to deal with the unusual non additivity condition, see the section devoted to equilibrium. 
The situation is somewhat similar as far as relaxation to equilibrium is concerned: in 
this case also, classical tools, namely those of kinetic theory, have proved sufficient to 
explain some unexpected phenomena. Thus, we fell that there is at present no obvious 



need for an alternative theory describing the relaxation of these systems with long range 
interacting. 

Let us note that despite the successes of these well established theories, standard sta- 
tistical mechanics at equilibrium, classical kinetic theory concerning the relaxation to 
equilibrium, there remains open questions and challenges, and room for new discover- 
ies, especially concerning the relaxation; we have tried to outline a few of them along 
the way. The most important are probably, on one hand the quest for ensemble inequiv- 
alence, negative specific heat and phase transitions in natural phenomena or laboratory 
experiment, and on the other hand the understanding of the limits of kinetic theories. 
However, we feel that the most relevant questions, both theoretically and practically, 
concern forced and dissipative systems, out of equilibrium. The last section presents 
very recent preliminary steps towards an undertanding of these situations, for which the 
theory is far less developed. 
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